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1 Introduction 

Forward-backward asymmetries in the production of fermions in high-energy e~^e~ coUisons 
are known to be precision obervables for the determination of the respective fermionic neutral 
current couphngs. Specifically the forward-backward asymmetry A^p^ of b quarks, which was 
measured at the Z resonance with an accuracy of 1.7 percent, led to a determination of the 
effective weak mixing angle sin^ ^VK,e// of the Standard Model with a relative precision of about 
1 per mille ^1 12] ^ notwithstanding the apparent discrepancy between this measurement and 
the determination of sin^ ^w,eff with similar precision from the left-right asymmetry measured 
by the SLD collaboration. (For a comprehensive overview, see P].) 

At a future linear e~^e~ collider '^, precision determinations of electroweak parameters will 
again involve forward-backward asymmetries. When such a collider will be operated at the 
Z peak, accuracies of about 0.1 percent may be reached for these observables J3J Ej. More- 
over, the top quark asymmetry A^^ will be experimentally accessible - a crucial tool for the 
determination of the hitherto unexplored neutral current couplings of this quark. 

The theoretical understanding of these observables, in particular those involving the heavy 
quarks Q = t,b, c, must eventually match these projected accuracies. The present theoretical 
description of Ap^ and Ap^ includes the fully massive next-to-leading order (NLO) electroweak 
13 E] and fully massive NLO QCD CHI IH] corrections. The next-to-next-to-leading order 
(NNLO) QCD corrections, i.e., the contributions of a^ to these asymmetries, were calculated 
so far only in the limit of massless quarks Q [T^ITHI IT^. To be precise, the forward-backward 
asymmetry of massless quarks is not computable in QCD perturbation theory, as was pointed 
out in J3]. It is affected in the limit rriQ — * by logarithmic final state divergences ~ In mQ 
resulting from the contribution of the QQQQ final state. (These terms are associated with the 
non-pertubative Q fragmentation function.) These logarithmically enhanced terms were taken 
into account in Ref. ^1], which is the most complete calculation of A^^p^ within QCD to date, 
done both with respect to the quark and the thrust axis. 

In view of the future perspectives for the b- and t-quark asymmetries at a linear collider, 
a computation of the order a^ contributions to Ap^ for massive quarks Q is clearly desirable. 
The NNLO QCD corrections involve three classes of contributions: (1) the two-loop corrections 
to the decay of a vector boson into a heavy quark-antiquark pair; (2) the one-loop corrected 
matrix elements for the decay of a vector boson into a heavy quark-antiquark pair plus a gluon; 
(3) the tree level matrix elements for the decay of a vector boson into four partons, at least two 
of which being the heavy quark-antiquark pair. 

In the limit of massless external quarks Q the QQ contributions to Ap^ vanish up to a 
non- universal correction of order a^ due to quark triangle diagrams |12^ . This is no longer the 
case for ttiq ^ 0. In this paper we determine class (1), i.e., the order a^ QQ contributions 
to ^^5, for arbitrary quark mass mq and center-of-mass (cm.) energy y^. For this purpose 
we use our recent results on the two-loop vector and axial vector vertex functions for massive 
quarks [2II1 Ell i22| which determine the amplitude of e^e~ ^ •y* , Z* ^ QQ to order a^ and 
to lowest order in the electroweak couplings. The contributions (1) from the two-parton final 
state, and (2) plus (3), i.e., those from the three- and four-parton final states, are separately 
infrared- finite, as will be discussed below. The latter can be obtained along the lines of the 
calculations of three-jet production involving heavy quarks ^3^1 El- The parity- violating 
part of da {3 jet) at order a^, which is a necessary ingredient here, was computed in jH 



However, a full computation of Ap^ ^""^ °" has not yet been done for massive quarks. 

The paper is organized as follows. In Section 2 we set up the formulae for determining Ap^, 
off and at the Z resonance, to order a^ in terms of the symmetric and antisymmetric cross 
sections for inclusive production of quarks Q. In particular we express the QQ contribution to 
ApQ by the one- and two-loop heavy-quark vector and axial vector form factors that determine 
the amplitude corresponding to the QQ final state. Moreover, we show that these contributions 
to the respective forward-backward asymmetry, which we denote by Ap^, are infrared-finite. In 
Section 3 we present the NNLO QQ cross section and the forward-backward asymmetry in the 
energy region a^ -C /5 -C 1, where P denotes the heavy-quark velocity. This result is applicable, 
for instance, to t quarks in the vicinity of the ti threshold. In Section 4 we evaluate Ap^ for 
b and c quarks at the Z resonance, and in the case of t quarks for c. m. energies above the ti 
threshold till 1 TeV. We conclude in Section 5. 

2 The Forward- Backward Asymmetry to Order a'^ 

In this paper we consider the production of a heavy quark-antiquark pair in e~^e~ collisions, 

e+e-^Yiq),Z*iq)^QQ + X, (1) 

where Q = c,b, t, to lowest order in the electroweak couplings and to second order in the QCD 
coupling as- To this order the following final states contribute to the cross section of inclusive 
QQ production, Eq. ((T)): the two-parton QQ state (at Born level, to order a^, and to order a^), 
the three-parton state QQg (to order a^ and to order a^) and the four-parton states QQgg, 
QQqq, and QQQQ (to order a^). 

The forward-backward asymmetry^ ApB for a heavy quark Q is commonly defined as the 
number of quarks Q observed in the forward hemisphere minus the number of quarks Q in the 
backward hemisphere, divided by the total number of observed quarks Q. The axis that defines 
the forward direction must be infrared- and coUinear-safe in order that ylpB is computable 
in perturbation theory. Common choices are the direction of flight of Q or the thrust axis 
direction. The forward-backward asymmetry can also be expressed in terms of the cross section 
for inclusive production of quarks Q. We have: 

with the antisymmetric and symmetric cross sections a a and as = cr defined by 

(^A= I -1 -dcos^- I rrfcost?, (3) 

Jq acosv y_i dcosw 

(rs= -^^dcosi3. (4) 

J__i dcosv 

Here "d is the angle between the incoming electron and the direction defining the forward 
hemisphere (in the e^e~ center-of-mass frame). When choosing the momentum direction of Q 

^We drop the superscript Q in the foUowing for ease of notation. 



or the thrust axis, the QQ contribution to Ape, which we compute in this paper, is of course 
the same. In the following d = Z(e~, Q). 

In analogy to its experimental measurement, Ape may be computed by determining the 
contributions from the final-state jets which, to order a^, are those of the two-, three-, and 
four-jet states. These contributions are separately infrared-finite, and Apg would not depend 
on the jet clustering algorithm employed when no phase space cuts are applied. Such a calcu- 
lation would require a jet calculus for massive quarks at NNLO in as which is, however, not 
available. (For massless quarks a NNLO subtraction method was recently developed |19J, cf. 
also references therein.) Yet, Ap^ belongs to the class of observables that can be computed at 
the level of unresolved partons. The two-parton and the three- plus four-parton contributions 
to the second-order forward-backward asymmetry are separately infrared (IR) finite, cf. fl^ ITi] 
and Section 2.2 below. This basic result will be exploited in the following. 

To order a^ the symmetric and antisymmetric cross sections receive the following contribu- 
tions from unresolved partons: 
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where the first number in the superscripts (i, j) denotes the number of partons in the 
final state and the second one the order of ««. Inserting (0) into © we get for A^^ 
order in a,: 
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Taylor expansion of © with respect to a^ leads to 
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where v4pB,o is the forward-backward asymmetry at Born level, and Ai and A2 are the 0{as 
and 0{a1) contributions normalized to Afb,o- 
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2.1 The QQ Contribution 

It is convenient to rewrite Eq. ((7j) as follows: 
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where the superscript {np) labels the number of partons n in the final state. Collecting the 
two-parton term from Eqs. (|H1), ©, and (fTUj) we get: 
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The remaining terms in Eqs. Q, ()10|) contribute to A^^^^ and A'^'^p\ 

As already stated above, both the two-parton and the three- plus four-parton contributions 
^(2p) g^j^d A^^P^ + A^^P^ respectively, are infrared (IR) finite, i.e., free of soft and coUinear 
singularites. We shall show this explicitly for A'^'^^^ at the end of this section. 

Next we express Api in terms of the VQQ vertex form factors {V = 'j, Z) which determine 
the amplitude of the reaction e~^e~ —>■ 7*, Z* -^ QQ to lowest order in the electroweak couplings 
and to any order in the QCD coupling; see Fig. 1. In this case the VQQ vertex Fq depends. 
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Figure 1: The amplitude e^e — » QQ in QCD. 
for on-shell external quarks, on four form factors: 
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[l^, 7^^]) "^Q denotes the on-shell mass of Q, and, for f = e,Q, 
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Here ej and T? denote the charge of / in units of the positron charge e and its weak isospin, 
respectively, and sw {cw) are the sine (cosine) of the weak mixing angle 9w The functions Fj 
and Gi denote renormalized form factors; the renormalization scheme will be specified below. 
Instead of using the Dirac form factor Fi we express in the following the (anti) symmetric cross 
section in terms of 



F^{s) = F^is) + F^is) . 
Neclecting the electron mass we find for the two-parton contributions to crA,s'- 
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Dz = [(■5— "^D^+^TiIrl] with nizi T z being the mass and width of the Z boson, (3 = 

the heavy quark velocity, and N^ = 3 the number of colors. Because we have put rrie = the 

form factor G2 does not contribute to Eqs. ((201) , (EH)- The last term in ^IU\\. which contains 

Fz, is of higher order in the electroweak couplings as compared with the first term. We will 

neglect that term in the following. 

Expanding the form factors in Eqs. ()20|) . ()2ip in powers of a^: 
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leads to the expansions of the two-parton contributions as exhibited in (jSj), i.e., 
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The contributions to cr^ can be further decomposed as follows: 
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where the superscripts 7, Z, and •yZ label the pure photon and Z contributions, and the 7 — Z 
interference terms, respectively. These terms are given by 
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where the convention Fa = ReF^ + inlmFa, Ga = ReGa + ivrlmGa (a = 1,2) is used; i.e., a 



factor TT is taken out of the imaginary part. The antisymmetric cross section is given by: 
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With these formuale the QQ contribution ()12|) to the forward-backward asymmetry is expressed 
in terms of the one- and two-loop form factors Eqs. (j^^ - (j23). 

The second order term in the expansion ()24p of the axial vector form factor, G\ , receives 
so-called type A and type B contributions. Type A contributions are those where the Z 
boson couples directly to the external quark Q [^, while the triangle diagram contributions, 
summed over the quark isodoublets of the three generations, are called type B [22] • (In the 
terminology of [14] these correspond to universal and non-universal corrections, respectively.) 



Among Eqs. ^ - ^ only a^g'''^' and a^'^ depend on G^' . With G^' = G\ 
we separate in these terms the type A and B contributions: 
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Thus ()14|) can be written: 
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If one evaluates, in the case of bb and cc final states, Eq. (j48|) exactly at the Z resonance and 
neglects the contributions from photon exchange, then: 
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to the forward-backward asymmetry is ultraviolet- and infrared- 



2.2 Infrared-Finiteness of Al^ 

The renormalized vector and axial vector form factors (|T7j) were computed in |2(H l2H l22] to 
order a^ within QCD with Nf massless and one massive quark Q, in a renormalization scheme, 
which is appropriate for the case at hand: the wavefunction and the mass of Q are defined 
in the on-shell scheme while as is defined in the MS scheme. The renormalized form factors 
still contain IR divergences which are regulated by dimensional regularization in D = 4 — 2e 
dimensions. At one loop, F^ and G\ contain 1/e poles due to soft virtual gluons. At 
two loops, the dominant IR singularities of F^ and Gl ' are of order 1/e^ due to soft 
and coUinear massless partons, while F2 has only 1/e poles. However, when inserting the 
renormalized form factors into the formulae for the order a^ and a^ contributions to Api the 
IR singularities cancel and we obtain a finite result. This is well-known for A[ ^ , c.f. jllj . 

This cancellation of infrared poles can be seen in a straightforward manner if one takes into 
account the universal structure of the infrared poles in the form factors, which can be expressed 
in terms of one-loop and two-loop infrared singularity operators for a massive quark- ant iquark 
pair. 
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This factorization of the infrared poles is a well-known feature of massless multi-loop am- 
plitudes, where it can be derived from exponentiation [221 121] ■ In massive QED, the same 
behaviour was observed long ago [25j. To the best of our knowledge, infrared factorization in 

massive QCD was established up to now only at the one-loop level [211; the above pole structure 

(1 2) 
of the form factors suggests that it holds at higher orders as well. Expressions for Jl,^ can be 

read off from the explicit pole structure of the form factors [201 1^ ■ Exploiting that 
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one finds immediately that the two-parton contribution (J14j) to the forward-backward asym- 
metry is infrared finite. However, it should be kept in mind that the two terms in (|T^ are not 
separately finite. 

Let us illustrate how this cancellation occurs in A2 by considering this expression at the 
Z resonance, i.e., by taking into account only Z exchange contributions and neglecting the 
photon and 7 — Z interference terms. When inserting the form factors into Eq. (J15p . all the 
leading 1/e^ singularities cancel, but a subleading divergence remains: 
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where y = {1 — /3)/(l + /3) and Cp = {N^ — l)/{2Nc). This subleading singularity proportional 
to Cp results from the real parts of the one-loop form factors, F{ and Gi , and from the 
contribution of a set of abelian-type two-loop diagrams (gluon ladder diagram and gluon vertex 
diagrams with quark self-energy insertions) to the real parts of the two-loop vector and type A 
axial vector form factors. The singularity ()52j) is removed by the second term in (fT^ . While 
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and thus the singular part of A2^i cancels the singularity in (jH2I)- This cancellation pattern 
remains also away from the Z resonance. 

3 Cross Section and A^g near Threshold 

In this section we analyze the cross section and Api for e~^e~ — i> 7*, Z* -^ QQ near the pro- 
duction threshold. Close to threshold, where the quark velocity /3 is small, the fixed order 
perturbative expansion in as of these and other quantities breaks down due to Coulomb sin- 
gularities, which must be resummed. In addition, large logarithms in f3 appear, which may be 
summed using the renormalization group [23 |2H1 applied in the framework of nonrelativistic 
effective field theory methods. However, in the energy region where a.^ -C /3 <^ 1, threshold 
expansions of observables, i.e., expansions in f3 to fixed order in as are expected to yield reliable 
results. In this region we derive compact formulae for the symmetric and antisymmetric QQ 
cross sections at NNLO. These expressions should be useful, especially in the case of top-quark 
pair production, for comparison of continuum results with results obtained directly at threshold. 

3.1 The QQ cross section 

We start with Eqs. (|^ - (|^ where the QQ cross section at NNLO, (Tnnlo = <^s 1 i^ 
expressed in terms of the 7, Z, and 7 — Z exchange contributions. Using these formulae we 



may write 



O'NNLO 



a 



(2,0,7) 



1 + 



a 



(2,0,Z) 



+ a 



(2,0,7Z) 



a 



(2,0,7) 



+ a?^-) 



1 + 



a 



(2,0,Z) (2,1,Z) (2,0,7Z) (2,1,7Z) 



0" 



+ a 



s 



al 



a 



(2,0,7)^(2,1,7) 



S 



a 



s 



(2,2,7) 



1 + — 



(2,0,Z) (2,2,Z) (2,0,7Z) (2,2,7Z) 



Denoting 

^(0,Ax) 

and recalling Eqs. (jHB|) . (jH7|) we get 



a 



s 



+ a 



s 



ai 



(2,0,7)^(2,2,7) 



a 



(2,0,Z) , ^(2,0,7Z) 



S 



+ a 



5 



a 



(2,0,7) 



(55) 



(56) 



O'NNLO 



a 



(2,0,7) J ^ ^ ^(0,Ax) ^ ^(2,1,7) 



1 + /\iOAx) 



+ 



a 



(2,0,Z) 



0" 



(2,0,7) VS 



(2,1,Z) _ ^(2,1,7) 



+.f^-'' 



1 _L ^(0,Ax) 



+ 



(J 



CT 



(2,0,Z) 

_/ (2,2,Z) _ ^(2,2,7) 

(2,0,7) V -^ •^ 



Putting 



a 



(2,0,Z) 



0" 



CT 



(2,1,Z) 



— a 



(2,1,7) 



the NNLO cross section reads: 

O'NNLO 



(2,0,7) ^ •^ 

5 

(2,0,Z) 

' (2,2,Z) _ ^(2,2,7) 

(2,0,7) V -^ -^ 



^^2,0,7) I ^^^(0,A.) +^^2,1,7) 



+ar-^ 



^^(g)^^(2,A.), 

^^(0,A.)j^^^(^^j^(Mx) 
_^(0,Ax.)j+^^(^^j'^(2,A.)| 



(57) 

(58) 
(59) 



(60) 



We now expand the expression in the curly brackets of Eq. (jUU]) up to and including terms of 
order /3°. Using the threshold expansions [2011^1221 of the one- and two-loop vector and axial 
vector form factors and putting the renormalization scale fi = uiq we obtain 

ONNLo = 4''°'"^{l + A(°'^-) + C^(g)A(i'^^)(l + A(°'^^'))+C^(g)A(i'^^) 

+C'^(^^J'a(2.^^)(1 + A(0'^^)) +C^(^)'a(2.^--)}. (61) 

The Born cross section a^' '^ is given in Eq. (jSHI)- We get for the terms A^^'^*^) and A^^'^^^ 
which arise from the one- and two-loop photon-exchange contributions, respectively: 

6C(2) 



/^(l,Ve) 



^ -8 + OiP) 



(62) 
(63) 
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with 



A 



(2,ye) 



+ 24(2(2) - 4C(2) (^_ - 8 In 2 + 4 ln/3j - 4C(3) + 



A 



(2,ye) 
TVA 



(2,ye) 
(2,ye) 



12(2(2) 48C(2) 

+0(/5), 
^(l^-Y^^(2^0+^^^2)(^-161n2-61n/5)-26C(3) 

+0(/5), 
lfl(21n(2;3)-^)-,f-,0(/5), 



39 



151 



32 



C(2) 



176 
~9~ 



0{[3') 



(64) 

(65) 
(66) 

(67) 



In the equations above, ({2) = 7r^/6, Ca = N^, Tji = 1/2, and Nj is the number of hght quarks, 
which we take to be massless. For the threshold expansions of the terms A*^*'"^"^-* (i = 0, 1,2), 
involving Z boson exchange, we find: 






ill W+(-f)' 

Dz 



[vivir 



^{4f ^^ 



z\'^ 



3-/3^ 



+ (-g) 



+ 2 1 



0(/5^) 



.(2,A.) 64C(2)m^^(a|)^ ^f + (gf )^] 

^ - (.^.2)^(4m^-m|)^ ^- + ^^^) 



5 y 




, (68) 

(69) 
(70) 



where vj, hf are defined in Eq. (|T8|) . While the multiplicative factor A'-'^'"^'^-' is given in exact 
form, the higher order terms are expanded to the appropriate order in (3 such that the NNLO 
cross section ()61|1 is obtained to next-to-next-to-leading logarithmic order (NNLL) in (3. We 
observe: 

(i) To NNLL in (3 the NNLO cross section is infrared-finite, i.e., the QQ cross section is equal to 
the total production cross section to this order. This is due to the fact that close to threshold, 
real gluon emission is suppressed by a velocity factor with respect to the QQ final state. The 
two- loop triangle diagrams, studied in ^^, do not contribute to the cross section to this order 
m(3. 

(ii) In the expansions of A*^^'^"^) and A^^'^'^^ IR divergences appear at 0{(3) and C(/?2), respec- 
tively. In the expansions of A*^^'"^^) and A*^^'"^^^ such divergences show up to order (3^ and /3^, 
respectively. 

(iii) The threshold cross section cr^NLO was calculated to NNLL in /3 in the above energy region, 
for pure vector exchange, previously in j^ (see also ISO])- Putting the axial vector couplings 
to zero in the above expressions and comparing with [21] we find agreement, 
(iv) The first and second-order Z-boson exchange terms A*^^'"^^^ and A*^^''^^'^ involving the axial 
vector coupling dg, are of order /3^ and /3°, respectively. A^^'"^^') is thus relevant for an analysis 
aiming at NNLO accuracy at threshold (counting a^ ~ /5). The analytic result for the axial- 
vector contribution at 0{af), which can be obtained to all orders in /3 from [211 122]; has not 
been given before in the literature. The term A*^^'"^^) is small compared to A^^'^^^ see Fig. |21 

Nevertheless, at a high luminosity linear collider with polarized e~ and e"*" beams one may 
eventually be able to disentangle the vector and axial-vector induced contributions to the ti 
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Figure 2: The terms A^^'^''^ (dotted line), A(°'^^') (dashed hne) and A^^'^"") (sohd hne), ex- 
panded to order /5°, as functions of P, for the case of ti production, using mj=172.7 GeV. 
/? = 0.2 (/5 = 0.5) corresponds to a center-of-mass energy of about 357 GeV (404 GeV). 



cross section. A numerical calculation of axial vector contributions in the context of Lippmann- 
Schwinger equations was performed in [HI], and the axial- vector contribution to the threshold 
cross section at next-to-next-to-leading logarithmic order was determined in PB| by a combi- 
nation of numerical and analytical calculations. 

For completeness we mention that a summation of all terms of the form a^//3" should be 
performed when the region of small heavy-quark velocities /3 ~ a^ is approached. The result of 
the resummation of these leading terms (and of part of the subleading terms) is the well-known 
Sommerfeld-Sakharov factor (see, for instance [22] )■ Subleading terms may be resummed by 
the renormalization group [2H] in the context of effective field theory methods. In this paper 
we are not concerned with these summation methods, as we consider only the region /5 ^ a^. 



3.2 Antisymmetric cross cection and A^j^ 

Next we perform the threshold expansion of the second-order antisymmetric cross section 



a 



(A) 
NNLO 



0"^^ , given in Eqs. ^I^ and PH|1 - PU)) . in the same manner as was done above for 



(^NNLO, using the results of j21 El 122 • We obtain to NNLL in (3: 



a 



{A) 
NNLO 



where o")^ is given in Eq. (jHH|) and 
A(Ai) - 6C(2) 

A(^'2) 



271 



A(^'1) + C^(^)'a(^'^)} 



(71) 



P 



- 6 + O (/?) 



C^Aif '^) + C^Ai,^^) + iV.T^Ai:^'^) + T^Ai,^'^) 



Al"'^^) 



(72) 
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with 



A<«. . l!^_!5|a),24m_4C(2,(f-fl„2.Hl„,)-fc(3)-.f 

+0(/3), (73) 

a!-^^ ^ lfl(f-^M2,)),4C(2)(f-fln2-41n,)-fc(3)-14 

+0 (/?) , (74) 

^(21n(2/3)-^)+4 + 0(/?), (75) 



Al 



(^,2) 



Aj.^'^^ = -fc(2) + f + 0(/3^). (76) 

The term Aj^ ' in Eq. ()72j) is the contribution of the triangle diagrams computed in |22] • It is 
infrared- and ultraviolet-finite. In the case of ti production it is given by 

aL^'^) = C(2) (l6 In 2 - I) - 8 In 2 + ^ In^ 2 + O {[3') . (77) 

Notice that the the second-order antisymmetric QQ cross section (j71|) is infrared-finite to NNLL 
in P and is equal to the total antisymmetric cross section in this order. The terms A^"^'^) and 
^(-4,2) i^gcome infrared- divergent to order P"^ and P, respectively. 

Finally, the second-order forward-backward asymmetry is given near threshold by 

Aps = Afb,oCfb, (78) 

where Cfb is the ratio of the curly bracket in Eqs. (ffT|) and of the curly bracket in (j?)T|) divided 
by (1 + A^^'"^"^')). To NNLL in /3 it is equal to the complete forward-backward asymmetry Ap^- 
In Fig. 121 we have plotted the forward-backward asymmetry Eq. (f7S|) to order a^ for ti 
production above threshold in the range 0.2 < P < 0.5, where Cfb is the ratio of two expressions 
expanded to NNLL in p. The top mass and the other parameters were chosen as given in 
Eq. ()79p below. A comparison with the exact second order asymmetry Ai^-L will be made in 
Section 4.2. 



4 Numerical Results 

In this section we compute the QQ contributions to Ape to order a^ using the results of Section 
2 and the analytic results for the one- and two-loop vector and type A and B axial vector form 
factors given in [201 El |22j. The numerical evaluation of the harmonic polylogarithms that 
appear in these expressions were made using the code of JHII. For b and c quarks the above 
formulae are evaluated at and in the vicinity of the Z resonance and for t quarks between 
2mt < a/s < 1 TeV. The quark masses, whose values we use are given below, are defined in the 
on-shell scheme while a^ is the QCD coupling in the MS scheme. When calculating ApB for b 
and c quarks around the Z resonance, as is defined with respect to the effective Nf = 5 flavor 
theory; that is, the top quark contribution to the gluon self-energy that enters the vector and 
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Figure 3: The forward-backward asymmetry Apg to NNLL above the tt threshold in the range 
0.2 < /3 < 0.5 for [i = rrit. 

type A axial vector form factors to order a^ is absent. The forward-backward asymmetry for 
top quarks is computed in the six-flavor theory with the corresponding six-flavor QCD coupling 
determined from the five-flavor coupling at the matching point fi = rrit. We use the following 
input values P: 



mc = 1.5 GeV, rrih = 5 GeV, rrit 
771^ = 91.1875 GeV, 
sin^ Ow = 0.23153 



= 172.7 ±2.9 GeV, 
Tz = 2.4952 GeV, 
0.1187. 



as [mz 



(79) 



The value of mass of the top quark is the recent CDF and DO average [32] . For h and c quarks 
the type B axial vector contributions are evaluated with the central value of rrit given in (f?^ . 
In the following we denote ApB evaluated to order a^ and a^, respectively, by: 



AfJ^Ka.) 
AfJ^\a's) 



A 
A 



{QQ) 

FB,0 

(QQ) 

FB.O 



A'; 



(QQ)\ 



j^{QQ) j^ j^{QQA) 



A 



{QQ,B) 



4.1 A^^ for b and c quarks at and in the vicinity of y^ = niz 

Let us first consider the b quark asymmetry. As it is to be computed for y/s ^ mz we can safely 
neglect the masses of the u,d, s, and c quarks which contribute to the second order form factors. 
(As already mentioned above, the t quark contribution to the gluon self-energy is decoupled.) 
The type B axial vector form factor G^ is non-zero due to the large mass splitting between 
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Abb) 

^FB,0 


^m 


^{bb,A) 


^{bb,B) 


4b^(«s) 


4b ^(«D 


/^ = 


mz 
2 


0.103128 


-0.000365 


-0.000084 


0.001147 


0.103090 


0.103200 


/^ = 


= mz 


0.103128 


-0.000326 


-0.000100 


0.000919 


0.103094 


0.103179 


^^ = 


-- 2mz 


0.103128 


-0.000295 


-0.000109 


0.000753 


0.103097 


0.103164 



Table 1: The bb contributions to A-p^ for bottom quarks at ^/s = mz- 



t and b quarks, and to very good approximation one may neglect in these triangle diagram 
contributions the mass of the b quark. Therefore we use 



,i2i,B), 



,(2I>,B) 



Gr'"\s) = Gr"^(s,mfe = 0,mfe = 0) - Gp ^(s,mt,mfe = 0) 



{2i,B), 



(80) 



when evaluating Eq. ()48|1 . respectivley ()49|1 . for the b quark. The functions on the right-hand 
side of (jHDI), whose second and third argument denotes the mass of the quark in the loop and 
the mass of the external quark, respectively, are given in [331 122] • (We use the notation of |23-) 
Putting /3 = 1 in Eq. (jl^ . our result for A),^' '{s = m|) agrees with that of [12]. 
Table ^ contains the values for the lowest order forward-backward asymmetry at the Z reso- 
nance, together with the bb contributions to first and second order in a^ for three choices of the 
renormalization scale /x. The photon and 7 — Z interference contributions, which are (on the 
Z resonance) of higher order in the electroweak couplings, are not taken into account. Tabled 
shows that the QCD corrections are dominated by the type B contributions: they are about 
three times as large as the order a^ and about nine times as large as the order a1 corrections. 
This is due to the fact that we are close to the chiral limit, as mf,/mz -C 1. In this limit the 
order a^ vector and axial vector form factors F^ , G\ , and the order a^ vector and type A 
axial vector form factors F{ , Gi ' become equal while the chirality-flipping form factors 



l(QQ) 



I iQQ,A) 



vanish. In this limit A^ and A2 vanish, too, as an inspection of the above formulae 

shows. Thus the QCD corrections are dominated by the type B contribution. As it turns out, 
it amounts to a correction of the lowest order asymmetry by only about one per mille. 



(bb) Abb,A) 



I {bb,B) 



In Fig. |3]the first and second order QCD corrections A[ , A2 ' , and A2 ' , evaluated 
for /i = mz, are shown between 88 GeV < ^/s < 95 GeV. Here the contributions from photon 
exchange are included. Again the QCD corrections are dominated by the type B triangle 
diagram contributions. Varying the renormalization scale in the range m^/2 < /i < 2mz 
changes these numbers only by a small amount. 

Next we consider the c quark asymmetry for ^/s ^ m,z- To very good approximation we 
can neglect the masses of the c and b quarks in their contribution to the gluon self-energy. 

(2i B) 

Here the type B axial vector form factor G\ is again determined by the large mass splitting 
between t and b quarks. In view of the convention adopted in Eq. (fTTj) . where the neutral 
current couplings of the external quark are factored out, we use now 



G\ 



{2E,B) 



Gi '{s,mt,mc 



0) 



Gi ^(s,mb 



0,m, = 0) 



(81) 



when applying Eq. (PH|) . respectively (jl^ . to the c quark. Again we put /? = 1 in these 
equations. Eq. (jHT|) is equal in magnitude but opposite in sign to Eq. (jHOJ. 
Table 121 contains the values for the lowest order c quark forward-backward asymmetry at the Z 
without the 7 and •j — Z contributions -, together with the cc contributions to first 



resonance 
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Ace) 
^FB,0 


Af^ 


^(cc,A) 


Af'^^ 


Aifias) 


Aifial) 


/i = ^ 


0.073592 


-0.000170 


-0.000060 


-0.002418 


0.073580 


0.073397 


l^ = mz 


0.073592 


-0.000152 


-0.000063 


-0.001938 


0.073581 


0.073434 


H = 2mz 


0.073592 


-0.000138 


-0.000064 


-0.001586 


0.073582 


0.073460 



Table 2: The cc contributions to ApB for charm quarks at ^/s = rriz- 

and second order in a^ for three choices of the renormahzation scale fi. The QCD corrections 
are dominated again by the type B term, which in this case is about two per mille of the leading 
order asymmetry. The increase by a factor of about two as compared to the b quark results 
from the fact that |ff | < \v^\, c.f. Eq. fH^ . 



In Fig. inithe first and second order QCD corrections A^ , A. 



(cc) . {cc,A) 



and A2 \ 



evaluated 



for yU = mz, are shown between 88 GeV < ^/s < 95 GeV, including the contributions from 
photon exchange. For c quarks, too, the QCD corrections are dominated by the type B triangle 
diagram contributions. 
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Figure 4: First and second order QCD corrections A\ (dashed), A2 ' (dotted), and A, 
(solid) for fi = mz in the vicinity of the Z resonance. 



,(bb,B) 



In Fig. El the bb and cc forward-backward asymmetries to order af, AU{al) and Ap^ (a^) 
are displayed between 88 GeV < ^/s < 95 GeV, using /i = mz- Figs. 010 show that the order 
af asymmetry is increased, in the case of b quarks, by about one per mille and decreased, in the 
case of c quarks, by about two per mille of the respective lowest order asymmetry in the whole 
energy range considered. Varying the renormahzation scale in the range mz/2 < yU < 2mz 
changes the asymmetries shown in Fig. El only by a very small amount. 

The complete forward-backward asymmetry to order a^ for b and c quarks at and in the 
vicinity of the Z resonance is dominated by the respective three- and four-parton contributions. 
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Figure 5: First and second order QCD corrections A^ (dashed), A2 (dotted), and A2 
(solid) for /i = mz in the vicinity of the Z resonance. 
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Figure 6: ApJ{a'^) (sohd) and Ap^ (a^) (dashed) for /i = mz in the vicinity of the Z resonance. 
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Figure 7: Leading order asymmetry A^^q for three values of the top quark mass. 



A complete computation of these terms to order a^ has not yet been done for nih, rn, 
Nevertheless, we expect that the respective results of 
limit, will not change dramatically. 



c 7^ 0. 

which were obtained in the massless 



4.2 ^pB for top quarks above threshold 

Finally we compute the ti contributions to the forward-backward asymmetry in the reaction 
e"'"e~ — > tiX, sufficiently far away from the pair production threshold in order that perturbation 
theory in ag is applicable. That is, the following results apply to events with t and i velocities 
/3 ^ as- As already stated above, ag is defined in the six fiavor QCD, with all quarks but 
the top quark taken to be massless. The value of as{n = mt) is determined from the input 
value ()79|). and values of a^ at other energy scales are obtained by two-loop renormalization 
group evolution. Most of the results below are presented for three values of the top quark mass: 
the present central and 1 s.d. upper and lower values 172.7 GeV, 175.6 GeV, and 169.8 GeV, 
respectively. The type B contribution to A^l, Eq. PHj) . is computed with the two-loop axial 
vector form factor 



G 



i2i,B) 



(s) = Gf'^'^\s,mt,mt) - Gf^'^\s,mb 



0,mt) 



(82) 



which is given in 

In the following we consider the energy range 360 GeV < y/s < 1 TeV. In Fig. [7| the leading 
order asymmetry A^Jq is shown for three values of the top quark mass. In Figs. |H1 and IHl the 

order a^ correction A\ is displayed for three values of the renormalization scale /i and fixed 
top quark mass, and for three values of rrit and fixed /i, respectively. The anologous cases are 
shown in Figs. ^1 and ^2 for the order a^ correction A2 ' . The triangle diagram contributions 
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Figure 8: Order as correction A[ for three values of the renormahzation scale fi, using mt 
172.7 GeV. 
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Figure 9: Order a^ correction A\ for three values of the top quark mass, using /i = y^. 
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Figure 10: Order a^ correction A2 ' for three values of the renormahzation scale /i, using 
mt = 172.7 GeV. 
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Figure 11: Order a^ correction A2 ' for three values of the top quark mass, using fi = ^/s. 
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Figure 12: Order a^ correction A2 ' for three values of the renormahzation scale n, using 
mt = 172.7 GeV. 



0.0002 



0.0001 



-0.0001 



-0.0002 



-i 1 1 r- 



-i r- 




m,= 169.8 GeV --- 
m,= 172.7 GeV — 
mt= 175.6 GeV -■ 



300 400 500 600 700 800 900 1000 

V5 [GeV] 



Figure 13: Order a^ correction A2 ' for three values of the top quark mass, using /i = ^/s. 
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Figure 14: Forward-backward asymmetry Apg (^s) for three values of the renormahzation scale 
/i, using rrit = 172.7 GeV. 
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Figure 15: Forward-backward asymmetry Ap-J{al) for three values of the top quark mass, usinj 
l^= y/s. 
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Figure 16: Forward-backward asymmetry to lowest, first and second order in as using rrit 



172.7 GeV and /i = v^. AfJ^ (dashed), A 
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A2 ' are given in Fig. ^Jand Fig. ^Jfor three values of fi and rrit, respectively. From these 
figures we conclude that the two-parton QCD corrections to the lowest order asymmetry are 



(tt) 



1 {tt,A) 



moderate to small for ^i > 400 GeV. At ^i = 400 GeV, A{ is about 3.3 percent while A. 
is about 2.4 percent. As expected, the relative importance of the order a^ corrections increases 
as the centre-of-mass energy approaches the threshold region: for y/s ~ 360 GeV, A2 is 
larger than A\ , signaling that perturbation theory in a^ is no longer applicable. Contrary to 
the case of b and c quarks at the Z resonance the two-loop type B contributions are two orders 
of magnitude smaller than A2 ■ 



Figs. El and El show the forward-backward asymmetry A 



(ti), 



a: 



for three values of the 

renormalization scale and three values of the top quark mass, respectively. The dependence of 
the second order asymmetry on fj. is small: changing /i from y/s/2 to 2^/s changes ApB (af), 
for fixed rrit, only by about 1 percent at ^/s > 360 GeV, and this dependence on /i decreases 
with increasing c. m. energy. 

In Fig. Elthe tt asymmetry is displayed to lowest, first, and second order in as- This figure 
shows that for c. m. energies sufficiently away from threshold the QCD corrections are under 
control. 

Finally, a comparison is made in Fig. El between the exact second order forward-backward 



asymmetry A 



(it), 



a: 



as given in Fig. Eland the values obtained from the near-threshold NNLL 



formula Eq. (f7H|) . For ^i < 360 GeV corresponding to /3 < 0.3 the deviation of the NNLL 
from the respective exact value is less than 5 percent. 
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Figure 17: The second order forward-backward asymmetry Ape (ttg) near threshold: exact 
values (dashed) as given in Fig. ^1 and the values obtained from the near-threshold formula 
Eq. dZHl) (sohd), using /x = m* = 172.7 GeV. 

5 Summary and Outlook 

Future experiments on heavy-quark production at a planned linear e^e^ collider aim at very 
precise measurements of the neutral-current couplings of these quarks. An important observable 
for this purposes is the forward-backward asymmetry J^^ in inclusive heavy quark production, 
e^e~ -^ 7*, Z* — > Q + X. The projected accuracies with which A^^ can be measured in future 
b or t quark production requires also the precise determination of these observables within the 
Standard Model. In particular, a computation of the order a^ QCD contributions to Ap^ for 
massive quarks is mandatory. 

In view of these perspectives we have calculated the contribution of the QQ final state to 
ApQ in NNLO QCD. As discussed above, and explicitly shown for the QQ final state, the 
contributions of the two-parton and of the three- plus four-parton states to the second-order 
forward-backward asymmetry are separately infrared-finite. We have provided formulae for 
the symmetric and antisymmetric QQ cross sections as and a a which yield Ap^ . These 
formulae hold for any center-of-mass energy. Specifically, in the energy region near threshold 
where the quark velocity (3 satisfies a^ <C /3 <C 1, we have expanded the order a^ and a^ QCD 
corrections to as and a^ to NNLL in (3. To this order in (3 the QQ cross sections are equal to 
the corresponding total cross sections. Therefore, an (analytic) expression is obtained for the 
forward-backward asymmetry Ap^ to order a^ and order /3 near threshold. 

Moreover, we have computed the two-parton forward-backward asymmetry Ap^ for b and 
c quarks on and near the Z-boson resonance and for t quarks for center-of-mass energies ^/s 
above threshold to 1 TeV. The two-parton asymmetry is determined by the heavy-quark vector 
and axial vector form factors. To order a^ the axial vector form factors receive besides type 
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A (universal) corrections also triangle diagram contributions resulting from the large mass 
splitting between t and b quarks. These triangle diagram terms dominate, for ^/s ~ mz, the 
QCD corrections from the QQ final state to A^p^ and Ap^. This is due to the fact that here 
one is close to the chiral limit. However, the complete order a^ QCD corrections to these 
asymmetries are dominated by the contributions from the three- and four-parton final states, 
which were calculated so far only for massless quarks ^3] . For top quarks the triangle diagram 
contributions to A^p^ are negligible compared to the type A corrections. These corrections 
from the ti final state to the lowest order aymmetry are moderate for large ^/s and increase in 
size towards threshold. The order a^ corrections are important, as the analysis in Section 4.2 
shows. 

We plan to determine in the near future also the contribution of the three- and four parton- 
final states to the order a^ forward-backward asymmetry for massive quarks. 
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